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The Diurnal Variation of Terrestrial Magnetism. 

By George W. Walker, A.R.C.Sc, M.A, F.R.S., formerly Fellow of 

Trinity College, Cambridge. 

(Received October 15, — Read November 13, 1913.) 

The diurnal variation of terrestrial magnetism has been the subject of 
experimental study for many years, and at a considerable number of 
observatories scattered all over the earth. But the co-ordination of the 
results, and the theoretical investigation of the physical cause or causes of 
the phenomena, have not made progress which seems to bear a reasonable 
proportion to the vast amount of observational data that has accumulated. 

As far as I know, Dr. Arthur Schuster's memoirs* constitute the most 
inspiring and systematic attempt to reduce this matter to scientific law and 
order. Although I have found it convenient to depart from the method 
pursued by Schuster, this change of procedure was in no small degree 
suggested by various remarks in his memoirs. 

The comparative lack of interest in the subject arises, I believe, in great 
measure from the difficulty (commented on by Schuster) of obtaining the 
experimental data in a form convenient for rational comparison. Observa- 
tories still continue to reduce their observations and to publish their results 
in a variety of ways, and, unless one is prepared to undertake a considerable 
amount of tedious arithmetical computation at the very outset, it is prac- 
tically impossible to obtain a comprehensive view of what the facts are. 

The Advisory Committee of Eskdalenmir Observatory recommended that 
the instruments should be arranged so as to record directly the geographical 
components of magnetic force, and, when the tabulated results for 1911 
began to take definite shape, it was, perhaps, natural that I should be 
interested in comparing the results with those of the only other observatory 
which at the time also recorded the geographical components, namely, 
Potsdam (Seddin). The comparison, which I made in the light of Schuster's 
conclusions, brought out points of such interest that I was induced to attempt 
to collect data from other observatories as to the deduced values for the 
geographical components. Following Schuster, the data I desired to obtain 
were the Fourier coefficients in the geographical components. I confined 
myself to the 24-hour and the 12-hour terms, because I share the view, 
entertained by many of those who realise the errors that arise - on the 

* ' Phil. Trans.,' 1899, vol. 180, and 1908, vol. 208. 
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experimental and observational side, that the shorter period terms are of 
very doubtful accuracy. 

I further confined my attention to the average yearly value, not because 
the seasonal variation is unimportant, but because I felt that there was a 
great danger of coming to grief in a multitude of facts. 

The data I have been able to obtain are collected in tabular form. If they 
appear somewhat meagre, I would point out that I had perforce to confine 
my attention to observatories whose results could, without very much 
arithmetical labour, be put in the form desired. I hope that the information 
disclosed by the table will encourage others to make a contribution to it, and 
thus to the progress of knowledge of the phenomena. 

In the table the quantities are those in the equivalent formulae 

c'i cos t + h\ sin t ~f a 2 cos 2 1 + h 2 sin 2 1 

an d ci sin (t -f u\ ) -f c 2 sin ( 2 1 + a^), 

where t is the local mean time at the station and the unit of magnitude is 
0*1 7 or 10~ 6 gauss. 

From information very kindly supplied by the directors of the various 
observatories, the following statement may be made : — 

1. Practically all days are used, not selected quiet days. 

2. The original curves are not smoothed, except in the case of the V results 
for Pola. 

3. The original curves gave records of II and D except at Eskdalemuir and 
Secldin, where N" and W were directly recorded. 

4. The hourly values are those at the exact hour, except at Seddin, where 
the hourly value is the estimated mean for an hour centering at the exact 
hour. 

I need perhaps hardly say that one could wish for data obtained in 
precisely the same way, and that before entering on a minute arithmetical 
computation, all data should refer to the same year. It was soon obvious 
from my inquiries that one would have to wait several years before such 
data could be obtained. 

Nevertheless, I think certain broad inferences may be drawn from the 
data collected in the table, and these seem to me of vital importance. 

It will be convenient to review briefly the main points in Schuster's 
memoirs. If we grant the main proposition, that it is correct to represent 
the variations by a potential function, the problem may be divided into 
two parts: — 

(1) The empirical determination of the potential function which represents 
the observations. 
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(2) The theoretical investigation of the physical causes that give the 
potential function so determined. 

These two divisions are substantially represented in Schuster's first and 
second memoirs. But the two divisions are not entirely independent, (1) 
must borrow from any hint suggested in (2), and (2) must be kept within 
the bounds indicated by (1). 

In the first memoir, Schuster dealt with observations from the four 
observatories — Pavlovsk, Greenwich, Lisbon and Bombay. Observing that 
in the west component the Fourier terms depended substantially on the 
local time (a feature which is in the main confirmed by our table), he 
showed how the potential could be calculated. Since the west component 
is of the form SO /sin 6Scb, where 6 is the co-latitude and <£ the longitude, 
it is only necessary to get the proper expression for W in Tesseral harmonics, 
and then the simple integration with respect to <f> gives, when multiplied by 
sin 6, the function O. 

His conclusion reached was that the principal parts of the mean diurnal 
variation for the year could be expressed by a potential function of the form 

■Q,fa — Ai sin 6 cos 6 cos (t + u{) -f A 2 sin 2 6 cos 6 cos (2t + a 2 ), 

where t is the local time or the equivalent Greenwich time plus the 
longitude, say (t f + </>). 

In the second memoir Schuster proposes to explain these terms as arising 
from electrical currents in the atmosphere, set up by the joint action of the 
permanent part of the earth's magnetic field and the mechanical oscillation 
of the atmosphere. Thus the barometric variation is associated with the 
magnetic diurnal variation. Certain formidable difficulties occur and are 
carefully noted by Schuster, but I need not comment on them. 

In attempting to deal with the seasonal variation Schuster supposes that 
the conductivity of the air depends on the sun's zenith distance, and, 
assuming a simple and tentative expression for the conductivity, investi- 
gates the potential function that arises. The analysis is necessarily 
exceedingly complicated. 

Let us now return to the empirical expression obtained by Schuster. 

A potential of the form 

n/a = Ai sin 6 cos cos (£' + <£ + «i) 

must give components to north 1ST and to west W 

N" = + Ai cos 2 6 cos (t' + <f> + ai), W = — A x cos sin (£' + $ + ai) ; 
and again the form 

ft/a = A 2 sin 2 # cos# cos(22'-f 2<£-fa 2 ) 
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gives components 

N" = +A 2 sin^(2~3sin 2 0)cos(2^ + 2(/> + a 2 X 

W = — A 2 sin20sin(2^ + 2(jf) + a 2 ). 

The following Table II facilitates the comparison of the formulae just 
obtained with the data collected in Table I : — 

Table II. 



Station. 



Sind. 



Cos0. 



Cos 20. 



Sin 20, 



Sin 0(2-3 sin 2 0). 



Pavlovsk 

Eskdalemuir ... 
Wilhelm shaven 
Potsdam ... 
PeBilt ... 

Pola 

Helwan . . . 
Bombay ... 
Batavia . . . 



0*505 
0*569 
0*594 
0*612 
0*614 
0*709 
0*867 
0*946 
*994 



0*863 
0*822 
0*804 
0*791 
0*789 
0*705 
0*498 
0*328 
•0*108 






•490 





•352 





•293 





•252 





•245 


-0 


•005 


-0 


•504 


-0 


•790 


-0 


•977 



0-871 
0-936 
0-956 
0-968 
0*969 
1-000 
0-864 
0-614 
•0-214 















-0 

-0 

-0 



•624 
•586 
•560 
•536 
•534 
•349 
•221 
•648 
•958 



If now we examine the data and consider W only, then, making allowance 
for the fact that the data are not all for the same year, we have, I think, 
strong support for Schuster's expression. Batavia is, however, abnormal. 

Again, taking the values for N only we have strong support for a formula 
of the type given by Schuster. 

But clearly for any one station the values of N" and W should give the 
same constants in the formula. This is not the case. The phase angles 
are not in agreement, and what is even more serious is that the observed 
amplitudes of the north component ci and c 2 are all too great as compared 
with the values computed from the potential function which represents the 
west component. The matter looks still more difficult if we compare the 
a h Si, a 2 , h 2 of N with the b\, a h h 2 and a 2 of W. 

I confess that the result was disconcerting, for one could see no flaw in 
Schuster's method, and the west values did, on the whole, depend mainly 
on local time. There was, of course, the possibility of a higher zonal 
harmonic, but this did not promise much help, and I was certainly unwilling 
to enter on much arithmetical computation. There remained the possi- 
bility of obtaining a function which would contribute to N" without 
contributing largely to W. This practically meant a term depending on 
time from some fixed meridian, as Schuster had indicated. 

After carefully studying Schuster's theory of the cause, with its ensuing 
complicated analysis, it appeared to me that one might reverse the whole 
process, that is to say, one may in a perfectly general way determine 
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potential functions which are differentially related to opposite sides of the 
earth. The process may stop as soon as we get the terms which represent 
the data, and then we may seek for the physical meaning of the terms. 
In mathematical terms we start with the primary potential 1/r. Any space 
derivative of this is also a potential. Thus a space derivative fixed with 
regard to the earth gives a contribution to the fixed part of the earth's field, 
while a space derivative fixed with regard to the sun gives a contribution to 
diurnal change. 

In this simple way expressions were formed and it was found that the 
terms sought for could arise. I think it will be convenient to postpone 
consideration of the manner in which they were obtained and to give first 
the results, 

For the 24-hour term the form of potential function at which I arrived 
may be written 

O = A sin 5 cos 5 cos (t + ot{) + B sin cos (t + #2) 

+ C (sin 2 — §) cos (£ + </>'- <£) --I- D sin 2 cos (t + <j> - <£'), 

wherein t is the local mean time (= t' + <j>) and 0, <f> are the co-latitude and 
longitude of the station, while <£' may be regarded as the longitude of some 
arbitrary meridian. 

Since the data do not all refer to the same year a minute numerical 
analysis would be out of place. Thus only round numbers were selected for 
the constants. By trial the following form was finally selected : — 

10 6 fl/a = (150 sin cos 5 — 50 sin 0) cos £ + ( — 70 sin cos 5-10 sin 0) sin t 

+ 65 (sin 2 5-f ) cos (t—<j>) + 5 sin 2 cos (t + <f>). 

This gives in units of 0*1 7 the following values for the north and west 
components : — 

N = (150 cos 20 -50 cos 5 + 70 sin 2 cos <j>) cos t 

+ (—70 cos 25—10 cos 5 + 60 sin 2 5sin</>) sin t, 

W = (- 70 cos 5- 10 -10 sin 5 sin <j>) cos t 

+ ( — 150 cos 5 + 50 — 10 sin 5 cos <p>) sin t. 

The numerical values computed are shown in Table III. 

Although the observed values are given in units of 0*1 7, it must be 
remembered that the accuracy in most cases does not exceed 1 7. 

I think it will be admitted that we have, on the whole, got a substantial 
representation of the data. The better agreement that might have been 
obtained for the European stations had, of course, to be sacrificed somewhat 
to <*et the curious features exhibited by the Bombay and Batavia data. 
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Table III. 



Station. 


N. 


W. 


V. (Provisional.) 


a x . 


K 


a i. 


K 


a v 


K 


Pavlovsk 


+ 83 
+ 77 
+ 70 
+ 64 
+ 65 
+ 32 
- 49 
-142 
-137 


-16 
-36 
-21 
-13 
-20 
+ 7 
+ 57 
+ 87 
+ 59 


-73 

-67 
-67 
-67 
-66 
-61 
-49 
-42 
-1.2 


-84 
-79 
-76 
-75 

-74 
-63 
-32 
- 2 
+ 69 


+ 18 
+ 20 
+ 24 
+ 26 
+ 25 
+ 32 
+ 29 
+ 4 
-74 


-50 
-35 
-42 
-45 
-42 
-45 
-38 
-21 
+ 13 


Eskdalemuir 


Wilhelmshaven 

Potsdam 


De Bilt 


Pola 


Bombay 


Batavia 



The potential function obtained must also account for the values of V, 
We are at liberty to suppose that any term in X2 arises from an external 
or interior source, and further that these contributions may differ in phase. 
We may, therefore, introduce a division of the constants that will give the 
correct value of X2 at the surface and at the same time account for the 
observed values of V. In fact, if we could depend on the values of V we 
have the important means of determining the external and internal propor- 
tions. But unfortunately the observations of V are poor and everyone who 
has really faced the experimental difficulties admits that the results are most 
unsatisfactory. In particular I know that the Eskdalemuir results for V are 
quite unreliable, and, judging from the other data, I think they must be 
regarded with some suspicion. The Potsdam and De Bilt results do, however, 
agree very well, and so I take them as a sort of standard. 

But clearly the position is a weak one, and therefore I do not feel justified 
in giving more than a general indication of what our formula would do. 

I therefore assume that the outside and inside contributions agree in phase, 
and as regards the second order spherical harmonics I adopt Schuster's result 
that the internal contribution is one quarter of the external. In the first 
order harmonic I assume that the contribution from internal source is nil. 

We thus get the formula for V. 

V= {150 sin 6 cos 6 -50 sin + (70 sin 2 0- 130/3) cos <j>} cos t 

+ { -70 sin 6 cos 0- 10 sin + (60 sin 2 (9-130/3) sin <£} sin t. 

The computed values, marked provisional, are entered in Table III. The 
results are, on the whole, in the right direction, and might be brought closer 
by taking a smaller proportion from an internal source, but I do not think it 
worth while to force the matter until better data are available. 

VOL. LXXXIX. — A. 2 H 
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We now consider the 12-hour terms. The form which I found it convenient 
to try was 

12 = A sin 2 # cos 8 cos(2£ + «i) + B sin 2 # cos(2£ + a 2 ) 

+ iCsin 3 ^sin(2^^-^--(/)0 + I)sin6>(|~sin 2 6>)sin(2^-<^ + (/) , ) 

+ E sin sin (2£-<£ + <£'). 

The first four terms involve harmonics of order 3 but the last is of order 1. 
I finally dropped the term in B, made E = — -fD and D = \G. It also 
seemed an advantage to make </>' = 30°. 
Thus the empirical form adopted was 

10 6 X2/a= {~60sin 2 6>cos^ + 40sin 3 6>sin(^>-30 o )}cos2^ 

-f 25 sin 2 8 cos 8 sin 2t, 
so that G = 60, D = 20, E = -16. 

The components to north and west thus become 

N = {-60 sin 8 (2-3 sin 2 6>) + 120 sin 2 8 cos 8 sin(<£-30°)} cos 2t 

+ 25 sin 8 (2 -3 sin 2 8) sin 2t, 
W = {6Osin20-8Osin 2 0sin(<£-3O o )}sin2£ 

+ {25sin26> + 40sin 2 (9cos(^-30 o )}cos2^ 
The computed values are shown in the following Table IV : — 



Table IV. 



fll J • 


N 




W. 


V. (Provisional.) 


Station. 
















A* 


h> 


a 2 . 


b 2 . 


a 2 . 


b 2 . 


Pavlovsk 


-37 


+ 16 


+ 32 


+ 52 


-34 


+ 8 


Eskdalemuir ......... 


-50 


+ 15 


+ 35 


+ 68 


-44 


+ 11 


Wilhelmshaven 


-47 


+ 1^ 


+ 37 


+ 68 


-45 


+ 9 ! 


Potsdam 


-43 


+ 13 


+ 39 


+ 67 


-47 


+ 10 


DeBilt 


-47 


+ 13 


+ 38 


+ 71 


-49 


+ 9 


Pola 


-32 


+ 9 


+ 44 


+ 71 


-58 


+ 10 


Helwan 


+ 12 


- 6 


+ 52 


+ 51 


-59 


+ 14 J 


Bombay 


+ 63 


-16 


+ 42 


-12 


-10 


+ 24 


Batavia 


+ 45 


-24 


+ 3 


-89 


+ 69 


+ 29 



Again, I think these numbers give substantial agreement with the data. 

Turning to the data for V, the numbers look more hopeful than in the case 
of the 24-hour terms, and the experimental errors that affect the 24-hour 
term are perhaps not so serious in the 12-hour term. Thus I was tempted to 
push the agreement somewhat further. 
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The data suggested that the main part for the northern stations arose from 
the term in sin 2 cos 6, but was chiefly in the cosine term cos 2t and practically 
nil for the sine term sin 2t. This can be met by assuming a phase difference 
between the outside and inside contributions. Thus I assumed 

10 6 a/a = A ^sin 2 6 cos0 cos(2£-ei) + B ~ sin 2 cos 6 cos (2*-e 2 ). 

a 6 r 4 

Hence we have 

A cos ei + B.cos ei = — 60, A sin ei + B sin e 2 = 25. 
To get the V results for Potsdam, we have the additional equations 
3A cos ei— 4B cos e 2 = —120, 3 A sin ei — 4B sin e 2 = 0. 

Hence A = — 53, B = — 14, ei = —15°, e 2 = —50°, which agrees with 
Schuster's result that; B is about 1/4 of A. 

Next the term in sin 3 # suggested the main, part of the terms at Batavia, 
and assuming for the moment that the harmonic was entirely of order 3, I 
took the form 

10 6 &/& = sin 3 sin (<£ - 30) |^ 

Hence A cos ei + B cos e 2 = 40, A sin ei + B sin e 2 = 0. 

The values at Batavia suggest taking 

3 A cos6i — 4B cose 2 = 80, 3 A sinei — 4B sin6 2 = 40. 
Hence A = 35, B = 8, ei = 9°, e 2 = -45°. 

This again makes B about one-fourth of A. 

But the harmonic was really made up of a main part of order 3 and a minor 
part of order 1, and the correction to V is thus found to be 

sin 6 { - 16 sec 45° sin (<j> + 15°) cos 2t + 16 sec 45° cos (</> + 15°) sin 2t} y 

where we have assumed that the first order term was entirely external. The 
net result is 

V = { -120 sin 2 6 cos (9 + 80 sin 3 6 sin (0-30°) -23 sin 6 sin (£ + 15°)} cos 2t 
+ {40 sin 3 6 sin (0-30°) + 23 sin 6 cos (0 + 15°)} sin 2t. 

The computed values in Table IY seem to give very satisfactory agreement with 
the data. 

I do not, of course, suggest that this empirical representation of the data 
in Table I is unique or final, but I do think it is clear that the Schuster terms 
alone will not co-ordinate the data, and the form we have obtained does go a 
very considerable way towards meeting the difficulties. 

2 h 2 
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We may now consider how the terms were arrived at, in such a way as to 

provide a clue to the physical meaning 
without binding oneself to any special 
theory. 

In the figure let Of, 0??, 0£, be rect- 
angular axes through the earth's centre 0, 
in the directions, earth's way, radius to the 
sun, and perpendicular to the plane of the 
ecliptic. Let OX, OY, OZ, be fixed rect- 
angular axes in the earth, but not rotating 
with it, OZ being the axis of rotation, 
and OX in the plane containing 0£ and OZ. co is thus the obliquity of 
the ecliptic and % an angle expressing the time of the year. 

The direction cosines of OX, OY, OZ, referred to Of, Orj, 0£, are 

OX =' cos (o cos ^, cos co sin ^, —sin &> ; 

OY = — sin^/, cos %> 0; 

OZ — sin ay cos %, sin w sin ^, cos co. 

Let 0\ <j>' be the polar co-ordinates of a point fixed on the earth, and let 0, <fc 
refer to any other point on the earth. 

Suppose we start with a potential XI — 1/r and differentiate along the 
direction 0', <j>\ We get a new potential function 



,a 



, 3 



Xlj = ( sin f cos <£>' ^- -f sin 0' sin 6' *- + cos 

\ ax T dy 



, 3\1 



r\ 



OZ/ T 



{sin f sin cos (<p — <£') + cos 0' cos 0}. 

7 

Since $ — <£' remains constant while the earth rotates, we simply have a con- 
tribution to the permanent field of the earth, representing a doublet with its 
axis in the direction 0\ <£'. Again, let us differentiate along 0?;. We get 

3 1 



Ho zzz 



drj r 



V 



1 

= — — {sin cos </> cos co sin x + sin sin <f> cos ^ + cos 5 sin o> sin %} 

— { __ (l— sin 2 co sin 2 y)^ sin cos £ + cos sin &> sin v}, 

where £ is the local time at the point of observation reckoned from midnight. 
The term in cos contributes nothing to diurnal variation, but contributes 
a seasonal term to the general field. It has a maximum in summer, a 
minimum in winter, and vanishes at the equinoxes. 
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The term in sin contributes to the diurnal variation. The amplitude has a 
semi-annual seasonal variation, rising to its maximum value 1 at the 
equinoxes, and falling to its minimum value cos co at the solstices. 

Averaged for the year it contributes a term sin cos t to the mean 
diurnal variation. 

Such a term is indicated by the data. We may, if we like, picture it as a 
uniform field in the direction of the sun, within which the earth rotates. 

In a similar way we may differentiate along Of, and so get a term 
sin sin t. 

Let us now carry the process of forming functions further by differentiating 
Hi along the direction 0?;. . We may write the result in the form 

fl 3 = -^{Ssin^' sinco sin^ sin# cos# cos(<£ — <£')— cos#' sinco sin^ (1 — 3 cos 2 0) 

— 3 cos f sin cos 0(1 — sin 2 co sin 2 %)$ cos (<£ + 1') 

— 3 sin 0' sin 2 # (1 — mi 2 co sin 2 %)$ cos ($-\-t f ) cos(^>— -<£') 
+ sin 0' (1 — sin 2 co sin 2 ^)* cos (<£'-N')}, 

where t' is Greenwich mean time, G' } <£' the co-ordinates of a point fixed with 
respect to Greenwich, and 0, 6 the co-ordinates of the station referred to 
Greenwich. 

The first two terms contribute nothing to diurnal variation, but only 
annual seasonal terms to the general field. The remaining terms contribute 
to diurnal variation, the amplitude having a semi-annual seasonal variation 
fluctuating between 1 and cos co. Averaged for the year we get terms in the 
mean diurnal variation. 

If 6' is we have a term sin cos 6 cos t, where t is the local time, 
Similarly, a term sin cos sin t arises by differentiating along Of. We 
thus get the Schuster terms in the 24-hour term of diurnal variation. 

If 0' = |-7r, we have the term 

3sin 2 0cos((£ + Ocos((£ — 0') — cos(£' + O- ( A ) 

Similarly, differentiating along Of, we get a term 

3 sin 2 sin (0 + 1') cos (cj> - <£') - sin (<£' + 1'). 

ISTow, in this latter form replace <// by |-7r-f-<£', as we are at liberty to do. 

We get 

3sin 2 ^sin(</> + ^X siri (^-^ / )~ C0S (^ + O- (B) 

Adding (A) and (B), we get 

(3sm 2 0-2)eos(0'-M')> 
where t f is Greenwich mean time. 
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Again, deduct (B) from (A) and we get 

3 sin 2 9 cos (2 <£-<£' + 0- 

These are the terms we found it desirable to introduce. 

We have thus obtained the terms required for the data, and clearly we 
might have taken a more general form had it been necessary. 

We have seen that Schuster proposed to explain his term as arising from 
the permanent field and tidal oscillation. But there are great difficulties in 
the way. Our analysis confirms the association with the permanent field, 
but otherwise may simply mean differential conductivity of the air as 
between midday and midnight, and as between sunrise and sunset, the 
difference arising from the ultra-violet radiation from the sun. 

As regards the term depending on Greenwich mean time, its mode of 
generation would appear to associate it with a magnetic axis lying in the 
plane of the equator. If there is. such an axis the intensity will have to be 
feeble if it is not to produce a very pronounced influence on the general field, 
as indicated by observations. 

My main objective has been the mean diurnal variation for the year, 
letting the seasonal variation take care of itself. But I anticipated that any 
correct method of approach would provide by extension for the seasonal 
change. 

This appears to be the case in the analysis given. Seasonal change is 
provided for, but the terms we have obtained contribute only semi-annual 
change, not annual change. 

The latter is important, but will not of course leave an influence on the 
mean yearly value. Thus it may be expected to arise in a way which 
differs from that we have already used. But differentiation along f supplies 
the kind of function wanted. 

Thus -tt-t; = * {cos co cos 6 — sin co sin sin(£—i|r)}, 

where t is the local time and tan i|r .= cos co tau >%. 

The first term contributes to the permanent field. 

The second term contributes to the diurnal variation at any particular 
time of the year. But the phase ^ passes through a cycle in the course of 
the year, and so contributes nothing when averaged for a year. 

The seasonal change is, however, beyond the limits I had prescribed 
myself in this paper, and I shall not pursue it farther on this occasion. 

So far the examination has referred to the 24-hour terms, but the 
data indicate the necessity for adding to the Schuster expression 
sin 2 (9 cos #cos(2£ + a 2 ) in the 12-hour terms also. Our method for the 
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24-hour term suggests that we should get new 12-hour terms by a second 
differentiation along Of or 0?;. The formulae become more complicated if 
we are to retain the sun's annual variation in declination. 

We therefore simplify the matter by regarding the sun as remaining in 
the plane of the equator. If we thus lose sight of the seasonal variation, we 
gain in conciseness. 

If we take the form tt-t, -, i& -?r^ > we get as a potential form %^-~ 

orj J r orj r 6 

which leads to 

-A {(1 -f sin 2 <9)-f sin 2 6 cos 2t}. 

The contribution to the 12-hour term is of the form sin 2 cos 2t, but I 
did not find this term of much assistance in representing the data. 

The Schuster term sin 2 6 cos 6 arises by putting Q' = in fli, which 
then becomes cos 6/r 2 , and then forming the d 2 /dr] 2 , 3 2 /9f 2 , or 3 2 /3f 3?7 of 
cos 0/r 2 . 

If, on the other hand, we put 6 r = \ir in fli, which then becomes 
1/r 2 . sin 6 cos (<£ — </>'), and then form the 3 2 /3?7 2 , 3 2 /3|; 2 , or 3 2 /3f3?7 of this, 
we arrive at 12-hour terms of the forms — 

3 sin 6 cos (<£ + <£' + 20— ir sin3 e cos 2t cos (<£-~<£') ( C ) 

and 3 sin 6 sin (</> + $' + 2t')— 1 -f- sin 3 6 sin 2t cos (<£-<£')• 

In the latter replace </>' by J 7r + </>'. We get 

3sin6>cos(^ + ^ / + 20-^ sin36,sm2 ^ sm (^ ) --^ / )- C°) 

Addition of (0) and (D) gives the form 

sin 6 (!— sin 2 6) cos (2 £ + </>'-</>), 

while, if we subtract (D) from (C), we get the form 

sin 3 cos (2 £ + <£ — <£'). 

I found it convenient to add to these forms a term sin 6 cos (2t + <£' — <f>) 
which may be regarded as a uniform field parallel to the plane of the equator 
rotating with the angular velocity of the earth but in the opposite direction. 

We have thus formed a simple specification of the terms which seem to fit 
with the data. I have no theory to propose as to their origin, beyond 
Schuster's view that they arise from differential conductivity of the upper 
regions of the atmosphere. If further data confirm the probable reality of the 
terms, it should not be a difficult matter to express the law of conductivity 
which would account for them. As we have only had to proceed to second 
differentials with regard to the sun's direction, or the earth's way, it does not 
look as if the differential conductivity required is of a very complex nature. 



39*2 The Diurnal Variation of Terrestrial Magnetism. 

The possibility of associating the electric and magnetic state of the earth 
with its translational movement has often attracted attention, although 
without success. It may not be out of place here to give a general solution 
of the electromagnetic equations for a body of conductivity c moving in a 
straight line with velocity JcC, where C is the velocity of radiation. The 
solution, which I believe is new, refers to the steady state. 

If (X, Y, Z) («, ft, 7) are the components of electric force and magnetic 
force the equations are 

/ ^±,VlV7\- (dy__dft da dy dft dx\ 
\ ox J \oy oz oz ox ox oyj 

-k~)(* 8 )=fil-^ ^-^ 3X-3X' N 
dx/ ' \ dz dy ' dx dy ' dy dx 

and ~-f-~ + ^=0, 

ox oy oz 

\ dx /\dx dy dzj 
A solution is expressed by 

X = 0, Y ==-*?&, Z = *?*2, 
where 1^0 is a solution of 



(l-^^ + fe^ + ^+Sfc^ 0. 

ox* ox Olf oz* 



Primary solution- 



1 hc x ± hc a 

where /o 2 = - — ^ + y 2 + 2 2 , 



1-F 
and other solutions may be obtained by the usual process of differentiation. 



